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FIGURE 1. Total field intensity over the state of Minas 
Gerais, Brazil. The contour interval is 100 
nT 


FIGURE 2. 


Total intensity anomaly over the state of Minas 
Gerais, Brazil. The contour interval is 10 nT. 



INTRODUCTION 


This investigation aims at understanding the deep structure 
of the continental crust in Brazil and its variation in different 
structural environments. It is expected that a preliminary 
crustal model for the main structural provinces will be obtained. 
A secondary objective of this investigation is to check the 
normal geomagnetic field for the Brazilian area. 

During the period covered by this report our activities 
were concentrated on the following topics: 

i) Development and test of software to convert Magsat data 
tapes to a form readable by our computer. 

ii) Assessment of the contents of the first Investigator B test 
tapes . 

iii) Definition of our future activities in view of a preliminary 
analysis of the first total intensity anomaly map received 
from GSFC. 

iv) Collection and organization of correlative data. 

In the following sections the topics listed above will be 
discussed in greater detail. 

DATA CONVERSION 


All data tapes we receive must be converted from the original 
IBM/360 binary format to the Burroughs/B-6700 binary format. 

This is done in a two step operation. First a general Algol 
program was written which contains the code necessary to handle 
the conversion of the IBM variables of types R*4, 1*4, 1*2, L*l, 
A*4 and also two peculiar forms of integer variables found on 



The Algol source program is the input to a FORTRAN program 
together with the description of the format of the records to be 
converted. The output of the FORTRAN program is an edited version 
of the Algol source tailored to the specific tape to be converted. 


This two step procedure produces a very efficient code and 
is also very flexible, permiting to handle both the Investigator 
B and Chronicle tapes, as well as any unexpected format | 

alterations, as for instance the Investigator B header record i 

change in effect since March S 1981. 

, ! 

1 

Prior to be put to use, the data conversion software was I j 

1 » 

subjected to intensive testing in order to identify «nd correct f ! ! 
possible programming errors. 1 


The listings of both the Algol source and the Fortran editor;, 

i ^ i 

are given in the Appendix A. 


DATA ASSESSMENT 


Two programs, LISTER1 and LISTER2, were written to provide 
formatted listings of the Investigator B tapes. 

Both programs read the (converted) Investigator B tape 
through the subroutine RDINVB. This routine performs all the 
I/O operations and also handles exceptional situations like 
end-of-file and I/O errors. The program LISTER1 produces an 
abridged listing of an Investigator B tape. From each header 
record only the Julian date, the pass number, the longitude of 
the nodes and the planetary activity index are printed. From 
each subsequent data record only the Julian date and time of the 
first data point are printed along with the pass number, the 

r i « o f n e f ho 1 toe A tt AY AV A7 ■ ♦> a«4 



This kind of listing is useful to verify the coverage of 
the data as well as the geographic selection and K index 

r 

screening of the data made by GSFC. 

The program LISTER2 produces an almost complete listing of 
an Investigator B tape. The header records are printed in full, 
except for the coefficients of the internal field model which 
are not printed. The contents of the data records are completely 
printed out with the adition of the anomalies AB, AX, AY, AZ 
which are computed, printed and printer-plotted. 

The listings of the two programs along with samples of the 
output are given in Appendix A. 

By examining the listings produced by the program LISTER2 
some inconsistencies were found on the first two INVB test tapes. 
At first it was suspected that the data conversion software was 
in error, but these inconsistencies were traced back to the 
original binary tapes received from GSFC. It was found that the 
components XMD and ZMD, computed from the spherical harmonic 
model coefficients. had the sign sistematically reversed with 
respect to the observed values for these components. Furthermore 
it was found that the field QUAL (attitude quality flag) specified 
as an 1*4 value in the description of the INVB format, was in 
fact an 1*2 value left justified in a 4 byte field. 


PRELIMINARY ANALYSIS OF THE GSFC GLOBAL ANOMALY MAPS 


A new version of an initial anomaly map, distributed to the 
investigators during the meeting at GSFC on December 4,5 1980, 
was received at the begining of May, 1981. 

As GSFC distributed only the map but not the gridded anomaly 
values used to plot it, we were restricted to a visual analysis 



The anomaly pattern in the region of South America is 
peculiar in two aspects. First, it is devoid of intense anomalies 
as compared for instance to the North American continent and to 
Australia. This is particularly visible in the colored version 
of the anomaly map. Secondly, the existing anomalies are weakly 
correlated with the main tectonic units of the Brazilian territory 
(DNPM, 1971b) except for the Amazon basin, which stands out as a 
prominent high with its maximum running along the basin. 

Consi ering that the potential resolution of the Magsat data 
is certainly better than 2°, as evidenced by the examination of 
the anomaly profiles of individual passes, we feel inclined to 
apply the equivalent source technique to the data. We hope that 
the cost in terms of computer time will be justified if the 
crustal magnetization parameters show such a high correlation 
with the tectonic pattern as was reported by Mayhew et al. (1980) 
for the POGO data. 

CORRELATIVE DATA 


A description of the correlative data we intend to use in 
the analysis of Magsat data, and which is adequatelly organized 
at present, is given below. 

i) Geological data 

The geological and tectonic maps of Brazil at the scale of 
1:5,000,000 are available (DNPM, 1971a, 1971b). 

ii) High resolution aeromagnetic data 

An aeromagnetic survey covering the region limited by the 
parallels 14°S and 23°S and the meridians 39°W and 38°W was 
digitized from maps at the scale of 1:100,000. The digitization 
was made on a regular 0,1° x 0,1° grid. 


Maps of the measured field intensity and of the total 
intensity anomaly, plotted from the digitized data, are shown 
respectively on figures 1 and 2. 

iii) High resolution gravity data 

A gravity survey covering an area delimited by the parallels 
18°S and 20°S and the meridians 41,5°W and 50°W is available. 

The Bouguer map of the western part of this survey, which 
overlaps the area covered by the aeromagnetic survey previously 
described, is given in the paper of Blitzkow et al. (1980). The 
Appendix B contains a copy of this paper. 

iv) Low resolution gravity data 

There are two sources of low resolution gravity information 
we intend to use. One of them is the gravity anomalies derived 
from satellite geopotential models like GEM-9 and GEM-10 (Lerch 
et al., 1977). The other source is the set of terrestrial 
measurements averaged over 5° x 5° blocks (Gaposhkin, 1979). 


PROBLEMS 


We have been experiencing difficulties with the customs I 

clearance of some of the postal parcels containing Magsat Data j 

•i 

tapes. 

We referred the problem to Mr. Locke Stuart, and he suggested 
that we should try to receive the Magsat data by diplomatic 
pouch through the Brazilian embassy in Washington D.C.. 

Unf ortunatel ly we cannot use the diplomatic service for this 
purpose. 

According to the customs authorithy, the problem is that 


involved, and this va ie is not clearly specified for the Magsat 
data parcels. In orde • to correct the problem, all parcels 
should be accompanied by either an invoice of zero monetary value 
or a statement issued by the sender declaring the contents of 
the parcel to be scientific material with no commercial value. 
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A gravity survey was carried out in 
the Slate of Minas Gerais, between 18° 
and 20°S, 42° and 44°30' W, utilising 
the FV3GE topographic network. The- 
se first order data were integrated with 
other stations where the altitude was 
measured by barometric methods. The 
x un'ey was made along profiles with the 


Foi efetuado um levantamenlo gra- 
vimetrico no Estado de Minas Gerais, 
entre '8°e20°S, 42°e44°33‘W. sobre 
a rede de nivelamenlo geomitrico da 
FIBGE. Nos inlervalos onde ndoforam 
encontradas KNs as altitudes foram 
determinadas por nivelamenlo baromi- 
trico. O levantamenlo gravimitrico foi 
amduzido ao longo de perfis e tern o 


ABSTRACT 

aim of extending the crustal informa- 
tions obtained by an unreversed refrac- 
tion seismic profile. Bouguer anomalies 
were inverted by computer techniques 
which make use of two-dimensional 
models of arbitrary shape and fitting 
with experimental data by iterative me- 
thods. 


RESUMO 

objetivo de eslender as informayoes so- 
bre a estrulura crustal, obtidas at r avis 
de uma prospecfdo sis mica a refrayao 
ndo revertida. A inversdo das anoma- 
lias de Bouguer foi efetuada com tec- 
nicas de computafio ulilizando mode- 
los bi-dimensionais de forma qualquer 
e ajuste uos dados experimentais por 
metodos iterativos. Foi obtida uma 


A normal depth of 38-40 km teas 
obtained for the Mono in the area. A 
crustal thickening was observed bet- 
ween Itabira and GuanhOes where the 
Moho reaches a depth of 45 km. . 

A lower crustal layer of intermediate 
density is required from gravity data. 


profundidade normal de 38-40 km para 
a MOHO. 

Entre Itabira e Guanhhes, foi obser- 
wda uma espessura crustal maior e a 
MOHO se encontra a uma profundidc- 
de de 45 km. 

Os dados gravimitricos implicam na 
exist fncia de uma camada com densi- 
dude intermedUtria na crosta inferior. 


Gravity data processed by modern 
techniques proved to be very useful la 
preliminary investigation* of tbe main 
crustal featurs and tbe crust-upper 
mantle transitions in different slruclu- 
ral environments. In many cases they 
furnished low cost preliminary crustal 
models, wbicb were afterwards 
i checked and detailed by more sophis- 
ticated and expensive techniques, such 
as deep seismic soundings (DSS). 
They ^ave been used also as a metho- 
dology complementary to l)SS to ex- 
tend to a larger area the crustal Infor* 


INTRODUCTION 


■nations obtained along single IWS 
profile*. 

Successful gravity derived c r u stal mo- 
dels exist in oceanic atom (e.g. Taiwan! 
et al., 19651, bland ai -s (e.g Makris, 
19761, young orogeas ( 04 . MorcW, 
1973* Corrode et tu„ 1974), pre-Cam- 
brtaa areas ( 04 . Ramberg, 1973). 


Quantitative estimates of the crustal 
thickness in Ike Brazilian territory are 
quite scanty at present. 

The aim of Ihb paper k to construct a 


gravity 1 

del la the isatbmlem part of the Mi- 
nas Gerais Slate (batsman 11* and 
J0*S, 43 and 44*30* Wl. ..mplamea 
ting and extending tbe laiarmattaaa 
f u rbi s he d by dm as wj DSS prods up la 
new existing hi that area. This b a 110 
km long onr o v en od prods which was 
carried out b e t ween I ta bi ra and Sle 


J 060 EvangeUsla by a German group 
directed by P. Gleet*, under the spon- 
sorship of the BrazU-Germany agree- 
ment for scientific cooperation (Giete, 
1975). 


GEOLOGICAL OUTLINE 

The surveyed tree includes the 
south-western margin of the Slo 
Francisco craton and the 
neighbouring metamorphic-sedimen- 
tary belt formed during the Brazilian 
cycle (Almeida, 1977), i.e. the latest 
orogenic cycle which affected a targe 
part of South American in late pre- 
Cambrian and Eo-Cambrian times 
(about 600-650 m.y.) (Cordani et al„ 
1973), 

A very simplified geological sketch 
map is reported in Fig. 1 High grade 
metamorphic terrains of granulite and 
amphibolite facies and kinzigitic 
gneisses are largely predominant. The- 
se formations are of Archean (2,000- 
2,600 m.y.) and pre-Cambrian (1,000- 
1,400 m.y.) ages and constitute the 
oldest outcropping rocks. A pre-Cam- 
')i»n metasedimentary cover (Espi- 
nhaco formation) overlies these rocks 
in a N-striking region which is located 
N of Quadrilatero Ferrifero. A N- 
trending reverse fault terminates the 
Espinhavo formation to the We$(, 
where it contacts cabonatitic rocks 
and metasediments belonging to the 
Brazilian cycle. 

The intermixing of rocks of diffe- 
rent metamorphic grade and chemical 
composition does not allow a reliable 
estimate of local densities to be used 
for Bouguer correction in the arias 
where the basement crops out. A 
slightly lower density might be used 
for the stations located on the Espi- 
nhaco formation and on Brazilian cy- 
cle terrains because in principle the 
sedimentary cover should have a 
lower density than the basement, ho- 
wev;r, as the sediments are generally 
metamorphosed and its thickness is 
unknown, it has been preferred to use 
the same density for the whole area in 
order to avoid contamination of the 
results by uncontrolled assumptions. 
Hence, the average density usually 
assumed for the upper crust (2.67 
g.cm*3) has been used for the Bou- 
guer corrections. No correlation bet- 
ween Bouguer anomalies and altitude 
has been found, even for stations lo- 
cated on the sedimentary cover. This 
indicates that the chosen density is a 
re sonable onr 


THE GRAVITY SURVEY 

About 3(X) gravity mcasuiemcnts 
Mere made in an area of 30,000 km2. 
Taking into account the t ransit abilit y 
of the region and the aim of the sur- 
vey, a distribution of the measure 
meni stations along profiles Mas pre- 
ferred to a regular netMoik A dtsiii 
bulion along profiles perpendicular to 
the mam geological estructures is parti 
cularly suitable lor inv crim^ .'a 


by computer techniques which make 
use of twodimeniional models. Three 
almost parallel NE trending profiles 
were chosen, as they area approxima- 
tely perpendicular to the main geolo- 
gical trends (from Curvclo to Couto 
de Magalhkes, from Belo Horizonte 
to SAo JoAo Evangelista, from Barlo 
de Cocais to Governador Valadares). 
Three transversal profiles (from Cur- 
veto to Sabari, from Datas to Gua- 
nhics, from SAo JoAo Evangelista to 
Governador Valadares) were also car- 
ried out. As the aim of the survey was 
not the investigation of shallow struc- 
tures a spacing of S km between the 
stations was considered adequate 

The gravity was measured by a 
Worden, Master type, and two La 
Coste and Romberg, Model G, gravi- 
ty meters. The measured values were 
referred to the base station of Pam- 
pulha Airport (Belo Horizonte), for 
which the IGSN value of: 

g - 978.385SO gal 

was used. Gravity values measured by 
Observatbrio Nacional (Rio de Janei- 
ro) (Gama, 1972) were also included 
in the construction of the Bouguer 
anomaly map (Fig. I) and the profiles 
(Fig. 3). These data were also reduced 
to the IGSN system by connecting 
some of the Observatbrio Nacional 
stations tc the Pampulha Airport base 
station. 

The gravity measurements were 
made on the FIBGE benchmarks 
available in the region. Further gravi- 
ty measurements were made on sta- 
tions where the altitude was determi- 
ned barometrically by System Pauling 
altimeters of Instituto de Geociincias 
of USP. by carrying out short circuits 
closed on FIBGE benchmarks. The 
accuracy of barometric determina- 
tions of altitude is estimated to be 
better than +■ 5m, leading therefore to 
an overall uncertainly not higher 
than 1 .0 mgal in the Bougher anoma- 
lies. 

Faye and Bouguer corrections were 
applied to the measured data. A den- 
sity of 2.67 g.cm'3 was used for the 
Bouguer correction (see discussion in 
the former section). Topographic cor- 
rections were not applied due to the 
lack of suitable topographic maps. 
Hovcwet station points were chosen 
generally in such a way as :o ensure 
that topographic effects were not 
higher than some tenths of a milligal 
as far as /one E of the Hayford 
scheme (360 m.). Topographic 
irrcgulaiities existing at longer distan- 
ces are not as prominent as to produ- 
ce el feet s higher than 0.2 mgal 

liouguer anomalies Merc computed 
Mith revpect to the 1V67 gravity refe- 
! e:.ce Meld. I he accuracy ol the Boo 
- icr anomaly tor each station is esti 


mated to be 1 .0 mgal, which is largely 
sufficient for the aim of the present 

ptpcf , 

The location of all gravity stations 
is reported in Fig. I, where open cir- 
cles indicate IBGE benchmarks, fitted 
circles indicate stations where the alti- 
tude was determined barometrically 
and irtenates indicate the stations 
measured by Observatbrio Nacional. 
A Bouguer map was elaborated (Fig, 
I). A listing of coordinates, observed 
gravity and reduction values for all 
the stations is available at Instituto 
Astronbmico e Geofisico, USP. 


THE BOUGUER ANOMALY 
MAP 

The prominent feature of the gravi- 
ty field in the area is a large low of 
about -50 mgal, which is located bet- 
ween Itabira and GuanhAes and is 
elongated in the same direction as the 
main geological structures (NNW). 
The extension of this anomaly indica- 
tes that its source must be very deep 
and it is likely that it is related to the 
crustupper mantle transition. The 
center of the .anomaly is located in the 
area where kinzingitic gneisses and 
amphibolitic facies rocks of the pre- 
Cambrian Barbacena formation are 
dominant. 

The pattern of the Bouguer anoma- 
lies in the surveyed area is quite regu- 
lar. The only evidences of important 
discontinuities in the upper crustal 
layers are a high of about + 30 mgal 
located SW of Serro and two lows, 
one located at the western margin of 
the surveyed area and the other bet- 
ween Governador Valadares and SAo 
Joio Evangelista. 

The steepening of the gravity gra- 
dient where the profiles C-C* and A- 
A‘ cross the contact between the 
EoCambrian SAo Francisco 
supergroup formations and the Espi- 
nhaco formations (Fig. 2) indicates 
the existence of a major tectonic dis- 
continuity. The occurrence of a simi- 
lar gradient around Itabira suggests 
that this tectonic line continues 
Southeastward, affecting also the area 
where the basement crops out. 

The data reported ,t. Fig. I were 
used for a quantitative evaluation of 
the depth of the Moho discontinuity 
ami any other possible intracrustal 
layer. 


CRUSTAL MODELS 

The Bouguer anomalies calculated 
along the profiles A-A\ B-B\ C-C" 
and D-D' (see Fig. I) were smoothed 
to eliminate the effects of sources lo- 


■EO-CAMSR'AN **0 FRANCISCO SURER OROUR (LIMESTONE, DOLOMITES, 
! SILTITCS, CALC0RHYLITES, AROCSIA*) 


RRE-CAMBRIAN COVER (ESRINHACO SUREROROUR QUARTZITES, RMTLUTES, 
M'TASEOWENTS, METAVOLCANIC) 



I ARE -CAMBRIAN AND ARCHCAR BASEMENT ( ORANULITES, AMRWBOLITES, GNEISSES) 


Ffg. 1 - faiyf m to trmtf mep rnd simplified g t otog toW tkath. C o n tam fcm knar val; i mgat. 


PROFILE A A 1 
• • PROFILE C C’ 


Y//A s FRSNCiSCO SUPERGROUP 
rxn ESPMHACO FORMATION 


*— -O 





catcd within ihe lew uppermost km of 
crust I tic obtained profiles were then 
inveticd by Ihe computer technique 
desetibed by Carrara et a! (1974). It 
consists ol the utilization of a two 
dimensional model ol arbitrary poli 
gonal cross section and assigned 
density contrast. At first,, the gravity 
field Dioduccd by the poligonal body 


the laylor expansion relative to the 
variable / (depth), limited to the 
terms of I s1 order. Then the compu- 
ted gravity is compared with the 
observed gravity by inicrativc pioce 
dures, which make use ol a least 
square method. The best fit i clause in 
the chosen model is so obtained. til 
ven ihe ambiguity inheient to giavity 


contrasts are not well known it is con* 
venient that at least the depth of the 
body be determined by other geophy- 
sical methods at some point along the 
prolile, if a reliable solution is 
desired. 

The only additional geophysical 
data existing in the area, which could 
help in defining the initial model lo 
be used in the computations, are (he 
P wave velocity functions obtained by 
the DSS unreversed profile (Gicse, 
197$). Two main P-wave groups were 
recorded along this profile. A wave- 
group with apparent velocities of 6.2 
km/sec was recorded as firts arrivals 
as far as the end of the profile (160 
k: i from the shotpoints). A second 
gioup had apparent velocity around 8 
km/sec. 1 tic long distance through 
which tlie first wave group is recorded 
as first arrivals indicates that the crust 
must have constant velocity dawn lo 
2S-30 km. No wave groups were ob- 
served which may indicate a first or- 
der disconti tuity within the crust. The 
analysis of the second wave group 
and the lacx of clear phases reflected 
from t lie Mt*ho discontinuity suggest 
that the t rut sinon (torn a velocity 
around 6.2 fm/sec to velocities typi- 
cal of ihe u(.pcr mantle must be con- 
tinuous and sharp Motto disconti- 
nuity exists. Y.»e depth where upper 






•I 42 km. I hcrclorc Ihe lowest 14 km • incraurt iik «|ipircni r-w«vc i iwu ainKnuww mmi wi inmiw; 
of crust must be characterized by P velocities shou*i be quite close to the shape indicates that: 

wave velocities gtadually incrcastni real velocities. Ihe crustal thickening a) the normal cruital thickness in 
with depth. occurs in a 80- 100 km wide NNW the surveyed area is 3I-4Q km; 


*M*IU *-«" 



The DSS profile coincides with the 
northermost part of the A- A' profile 
(from Cuanhles to Sko J oko Evange- 
lista). The crustal model used to 
invert gravity data is represented in 
Fig. 3. It was constrained to a crustal 
thickness of 42 km at Guanhkes. The 
chosen densities are consistent with 
similar models in other pre-Cambrian 
areas (e.g. Ramberg, 1973) and with 
the densities which can be estimated 
from apparent P-wave velocities 
applying the well-known Nafe and 
Drake relationship. Of course the 
assumption of a well defined lower 
crustal layer is a schemaii/aiion of the 
gradually changing layer suggested by 
DSS data. The density of 3.0 g.cnr’ 
is consistent with the average velocity 
of 7.4 km/sec suggested by Ciiesc 
(1975). The successful models and the 
computed profiles are reported in Fig. 
3a, and b. Standard deviations bet- 
ween computed and observed values 
are lower than ± 2.0 mgal for all the 

profiles. It must pointed or that 
the observed data are not compatible 
with any reliable two-layer-. model, 
that is the csisicncc of a lower crustal 
layer oi intermediate density i*. requi- 
red to have a good lit with the obser- 
ved data. Ihe satisfactory fit of the 
computed with the observed gravity 
indicates that the lov existing between 
Itahtia ami liii.mkVs is to be as- 
ciihcd pail Is to a deepening of Ihe 
Molio and put i Is to a deepening ol 
the top of the lowei crustal layer. I he 
crust icaches a iitastmuiii thickness ol 
43 km til ih;s aica. wheicas ihe nor 
mal depth all .inniiiil r- 3h 40 km 
I lie inciin.ii u>m o| ill. M.ilio along the 
seismic ptolilc i- mu luelici ili.m 3 


elongated area and it must be related 
with the compresxional tectonic boun- 
dary between the Focambrian Sko 
Francisco supergroup ar.d the 
Espinhaco formations. The highest 
crustal thickening occurs where me- 
dium and high grade metamorphic 
rocks of the Barbaccna formation 
outcrop. The crustal thickening repre- 
sentsprobably the remnant of an 
isoslalically adjusted pre-Cambrian 
mountain root. The constant P-wave 
velocity and the gravity derived model 
consistently indicate constant physical 
characteristics of the crust down to 
23-30 km of depth. No major varia- 
tions of metamorphic grade and/or 
chemical composition are therefore li- 
kely in this depth interval. The nor- 
mal thickness and the general charac- 
teristics of the crust in the surveyed 
area are close to those observed in 
other pre-Cambrian areas, such as 
Canada (Wynne- Edwards, 1972) and 
Southern Norway (Kancstr/tm, 1977). 

The lower crustal layer required 
from the gravity data begins from a 
depth of 20 to 27 km, according to 
ihe schematized model. This layer 
may be formed by recks of mafic 
nature. This inference is supported by 
the analysis of the long wavelength 
crustal anomalies of the magnetic 
field which suggest the existence of a 
compositional variation toward mafic 
terms at depth of 22 - 27 km. (Cor- 
rado cl a!., in press). 


CONCLUSIONS 

tin imi'tottu ol llicgiustls data by 
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b) a tone of crustal thickening 
exists between Itabira and GnanMtes. 
It is a 80-100 km wide NNW trending 
strip where the crust reaches • thick- 
ness of 43 km. The occurrence of a 
thicker crust in this area is consistent 
with the -xistence of a compressional 
front between the Sio Francisco 
supergroup and the Espinhaco forma- 
tions. 

c) the existence of a lower crust of 
intemediate density is required from 
gravity data. DSS data suggest howe- 
ver that no first order discontinuity 
exists within the crust, and that the 
italic characteristics of Ihe rocks 
constituting the lowermost crust 
approach gradually upper mantle 
values. If this situation can be sche- 
matized by assuming the existence of 
a lower crustal layer with density of 
3.0 g/cirw. it-.- depth of the lop of 
this layer is normally at 20 km, in- 
creasing to 27 km in the zone of crus- 
tal thickening. Greater depths can he 
computed for the top of this layer if 
higher densities are assumed. 
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